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Abstract 

Given a right-angled triangle of squares in a grid whose horizontal 
and vertical sides are n squares long, let N(n) denote the maximum 
number of dots that can be placed into the cells of the triangle such 
that each row, each column, and each diagonal parallel to the long 
side of the triangle contains at most one dot. It has been proven 
that Nf(n) = [ 2 " 3 ^ 1 j • In this note, we give a new proof of the upper 
bound Nf(n) < [ 2 " 3 ^ 1 j using linear programming techniques. 



1 The problem 

Consider a 'triangle' of squares in a grid whose sides are n squares long, as 
illustrated by the diagram in FigureQ] for which n = 7. We call a southwest- 
to-northeast diagonal a standard diagonal. Note that our triangle consists 
of all the cells in an n x n square that lie on or below the longest standard 
diagonal. 

We denote by N(n) the maximum number of dots that can be placed into 
the cells of the triangle such that each row, each column, and each standard 
diagonal contains at most one dot. Determining N(n) is equivalent to 
solving the following problem: Suppose we have a chessboard made up 
of hexagonal cells arranged in the shape of an equilateral triangle of side 
n. Then N(n) is the maximum number of non-attacking queens that can 



Figure 1 : A triangle of size 7 



be placed on such a board, where a queen can move in any one of the 
three directions allowed on a hexagonal grid. The following theorem was 
proven by Vaderlind, Guy and Larson [21 Problem 252] and independently 
by Nivasch and Lev [2]: 

Theorem 1.1. N(n) — Nf(n), where 
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Note that the value of Nf(n) can be stated more succinctly as follows: 
N f (n) 

In order to prove Theorem 11.11 we require a construction to establish 
the lower bound N(n) > Nf(n) as well as a proof of the upper bound 
N(n) < Nf(n). In [H [3], the upper bound was proven by elementary 
combinatorial arguments. In this note, we give a new proof of the upper 
bound using linear programming techniques. In the end, our proof is also 
combinatorial; the main contribution we make is to demonstrate the use of 
linear programming techniques in deriving the proof. 



2 The lower bound 



Before proving the upper bound, we give a construction to show that 
N(n) > Nf(n). This construction is essentially the same as the ones in 



11121- 
Theorem 2.1. N(n) > N f (n). 



Proof. First, we show that N(3t + 1) > It + 1: 



1. Place a dot in the leftmost cell of the (2t+ l)st row (where we number 
rows from top to bottom). 

2. Place t more dots, each two squares to the right and one square up 
from the previous dot. 

3. Place a dot in the (t + 2)nd cell from the left in the bottom row. 

4. Place t — 1 more dots, each two squares to the right and one square 
up from the previous dot. 

It is easily verified that at most one dot is contained in each row, column, 
or standard diagonal. 

Next, N(3t + 2) > N(3t + 1) > 2t + 1 (it suffices to add a row of empty 
cells). Finally, N(3t) > N(3t + 1) - 1 > 2t (delete the bottom row of cells 
from a triangle of side 3t + 1, noting that any row contains at most one 
dot). □ 

Example 2.1. We show in Figure [2] that A^(7) > 5 by applying the con- 
struction given in Theorem 12. II 

3 A new proof of the upper bound 

The computation of N(n) can be formulated as an integer linear program. 
Suppose we number the cells as indicated in the Figure |3] (where n = 6): 

Define Xi j = 1 if the corresponding cell contains a dot; define Xij = 
otherwise. The sum of the variables in each row, column, and standard 
diagonal is at most 1. This leads to constraints of the form 
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Figure 2: AT (7) > 5 



Figure 3: Labelling the cells in a triangle of size 6 
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and 

n 

^ £j,i-fe < 1, for k = 0, 1, . . . , n — 1. 

i=fc+l 

Finally, Xij G {0,1} for all The objective function is to maximize 
^2xij subject to the above constraints. It is obvious that the optimal 
solution to this integer program is N(n). 

It is possible to relax the integer program to obtain a linear program, 
replacing the condition Xij £ {0, 1} by < Xij < 1 for all In fact, we 
do not have to specify Xij < 1 as an explicit constraint since it is already 
implied by the other constraints; it suffices to require < Xij for all 
Denoting the optimal solution to this linear program by LP(n), we have 



Figure 4: The optimal solution to LP(6) 
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that LP(n) > N(n). 

Example 3.1. Using Maple, it can be seen that LP(6) = 4| (a solution 
to the LP having this value is presented in Figure [4} . 

Next, we tabulate some solutions to LP(n) for small values of n in 
Table [T] Based on the numerical data in Table [TJ it is natural to formulate 
a conjecture about LP(n): 

Conjecture 3.1 (LP Conjecture). Define 

LP f (3t) = 
LP f (3t + l) = 
LP f (3t + 2) = 

Then we conjecture that LP(n) = LPf(n). 

It is easy to show the following: 
Theorem 3.2. If the LP Conjecture is true, then N(n) = Nf(n). 
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Proof. First, the LP Conjecture asserts that 

LP(n) = LP f (n). 



(1) 



Table 1: Optimal solutions to the integer and linear programs for small n 
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Because iV(n) is an integer and N(n) < LP(n), we have that 

N(n) < [LP{n)\. (2) 
Simple arithmetic establishes that 

[LP f (n)\ = N f (n). (3) 
Combining (JTJ, @ and (|5J|, we have 

N(n) < [LP(n)\ = [LP f (n)\ = N f (n). 
We showed in Theorem HO that N(n) > N f (n); hence N(n) = N f (n). □ 

The optimal solution to the linear program for n = 6 that we presented 
in Figure |4] does not seem to have much apparent structure that could 
be the basis of a mathematical proof. Indeed, most of the small optimal 
solutions that we obtained are quite irregular, which suggests that proving 
the LP conjecture could be difficult. We circumvent this problem by instead 
studying the dual LP and appealing to weak duality. 

An LP in standard form is specified as: 



maximize c T x 

subject to Ax <b,x>0. 



This is often called the primal LP. Any vector x such that Ax < b, x > 
is called a feasible solution. The objective function is the value to be 
maximized, namely, c T x. 

The corresponding dual LP is specified as: 



Here, a feasible solution is any vector y such that A T y > c, y > 0. The 
objective function is b T y. 

We will use the following classic theorem. 

Theorem 3.3 (Weak Duality Theorem). The objective function value of 
the dual LP at any feasible solution is always greater than or equal to the 
objective function value of the primal LP at any feasible solution. 

We now describe the dual LP for our problem. Suppose we label the 
rows of our triangle by n, . . . , r n , such that is the row containing i 
squares, and we label the columns and diagonals similarly. The following 
simple lemma is very useful. 

Lemma 3.4. If a cell is in row r i7 column Cj and diagonal dk, then i +j + 
k = 2n + l. 

In fact, it is not hard to see that there is a bijection from the set of 
n(n + l)/2 cells to the set of triples 



In the dual LP, the variables are n, r 2 , . . . , r„, ci, C2, . . . , c n , d±, c?2, ■ ■ - ,d r 
There is a constraint for each cell C. If C is in row column Cj and diag- 
onal dk, then the corresponding constraint is 



The objective function is to minimize ^ f"j + ^2 c j + S ^fe • 

It turns out that there exist optimal solutions for the dual LP that have 
a very simple, regular structure. These were obtained by extrapolating 
solutions for small cases found by Maple. 

When n = 3t + 1, define 



minimize b T y 

subject to A T y > c, y > 0. 



T = {(i,j,k) : i+ j + k = 2n+ 1,1 < i,j,k < n}. 



n + Cj + d k > l. 





When n = 3t + 2, define 



n = c t = d l = max \ 0, -^—^ \ ■ ! : > 1 



When n = 3t, define 



d = di = max \ 0, ^ — ^ . (<> i 



Lemma 3.5. The values Ti,Ci anddi defined in (JS|) and © are feasible 
for the dual LP, and the value of the objective function for the dual LP at 
these solutions is LPf(n). 

Proof. First we consider the case n = 3t, + 1. Consider any cell C, and 
suppose C is in row r^, column Cj and diagonal d%. Then we have that 

i-t-1 j-t-1 k-t 

r, + ci + du > 
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(applying Lemma 13. 4 



Therefore all constraints are satisfied. The value of the objective function 
is 
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= LP/(3f + l). 

The proofs for n = 3t + 2 and n = 3t are very similar. □ 



Our new proof of Theorem 11.11 follows immediately from Lemma l3.5l bv 
slightly modifying the proof of Theorem 13.21 



Proof. First, from weak duality and Lemma 13.51 we have 

LP{n) < LP f (n). 
Combining this inequality with and ([3]), we have 

N(n) < [LP(n)\ < [LP f (n)\ = N f (n). 
We showed in Theorem |2~T1 that N(n) > Nf(n); hence N(n) = N f (n). □ 

4 Discussion 

We investigated the "dots in triangles problem" due to an application to 
honeycomb arrays (see [1]). However, we did not realize that the dots in 
triangles problem had already been solved. Since we did not know the value 
of N(n), we adopted an "experimental" approach: 

1 . We used Maple to gather some numerical data. 

2. We formulated (obvious) conjectures based on the numerical data. 

3. Finally, we proved the conjectures mathematically. 

Many problems in combinatorics are amenable to such an approach, but 
this particular problem serves as an ideal illustration of the usefulness of 
this methodology. Indeed, the problem seemed almost to "solve itself", 
with minimal thought or human ingenuity required! 

It should also be emphasized that, in the end, the resulting proof is 
quite short and simple: 

1. By a suitable direct construction, prove that N(n) > [ 2 " 3 ^ 1 j ■ 

2. Show that the dual LP has a feasible solution whose objective function 
value is less than [^^"J + 1- 

The first conjecture we posed was the LP Conjecture, concerning the 
optimal solutions to the LP. In general, to prove a feasible solution to an 
LP is optimal, it is necessary to do the following: 

1. Find a feasible solution to the primal LP and denote the value of the 
objective function by C. 



2. Find a feasible solution to the dual LP and denote the value of the 
objective function by C* . 

If C = C* , then the solution to the LP is optimal (this is often called strong 
duality) . 

When n = 1 mod 3, our work in fact proves the LP conjecture. This 
is because Theorem 12.11 yields a solution to the primal LP whose objective 
function value matches the solution we later found to the dual LP. How- 
ever, when n ^ 1 mod 3, we do not have a general solution to the primal 
LP whose objective function value matches the solution to the dual LP. Al- 
though we are confident that the LP conjecture is also true for these values 
of n, proving it could get messy! 
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